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Abstract 

The effect pointed out by A. B. Migdal in the 40's (hereafter named Migdal 
effect) has so far been usually neglected in the direct searches for WIMP Dark 
Matter candidates. This effect consists in the ionization and the excitation of 
bound atomic electrons induced by the recoiling atomic nucleus. In the present 
paper the related theoretical arguments are developed and some consequences of 
the proper accounting for this effect are discussed by some examples of practical 
interest. 
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1 Introduction 

The Migdal effect is known since long time and is described both in devoted papers 
[1] and in textbooks [2]; it has also been recently addressed for the Dark Matter 
(DM) field in ref. [3]. This effect consists in the ionization and the excitation of 
bound atomic electrons induced by the presence of a recoiling atomic nucleus. In 
this paper, it will be accounted in the case of the WIMP-nucleus elastic scattering. 
In fact, since the recoiling nucleus can "shake off' some of the atomic electrons, an 
electromagnetic contribution is present together with a recoil signal in the analysis 
of DM direct searches (with whatever approach) when interpreted in terms of WIMP 
candidates. Since this contribution is not quenched, one can expect that this part 
(usually unaccounted) can play a role as well. 
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In this paper the related theoretical framework is developed and some of the pre- 
vious corollary analyses for WIMP candidates [4, 5, 6] from the DAMA/Nal annual 
modulation results (total exposure of 107731 kg x day) are used for template purpose. 
We just remind that, in order to investigate in a model independent way the presence 
of DM particle component(s) in the galactic halo, DAMA/Nal [7, 8, 4, 5] has exploited 
over seven annual cycles the DM annual modulation signature, achieving a 6.3 a C.L. 
model independent evidence [4, 5]. Some of the many possible corollary quests for the 
candidate particlc(s) have also been carried out so far mainly focusing various possi- 
bilities for the class of DM candidate particles named WIMPs and for the class of light 
bosons [4, 5, 9, 6]; other corollary quests are also available in literature, such as e.g. in 
refs. [10, 11, 12, 13, 14, 15, 16], and many other scenarios can be considered as well. 



2 Atomic effects due to nuclear recoils 

The possible excitation and ionization of the atom by a recoiling nucleus - induced by 
a WIMP-nucleus elastic scattering - give rise to a certain quantity of electromagnetic 
radiation made of the escaping electron and of X-rays and/or Auger electrons arising 
from the rearrangement of the atomic shells. This radiation is fully contained in a 
detector of suitable size. 

Since the WIMP-nucleus interaction is expected to have a very short range (e.g. 
being mediated by very heavy particles with mass M), the duration of the collision (~ 
h/M < 10~ 26 s) is negligible with respect to the electron orbit periods and to R a /VA, 
where R a is the atomic size and Va is the nucleus velocity after the interaction. Thus, 
the interaction can be considered as instantaneous and a semiclassical description of 
the recoil process can be applied following the Migdal approach [1, 2]: the target 
nucleus is assumed to be at rest for t < t* and at t = t* it suddenly acquires the Va 
velocity because of the interaction. Therefore, after the collision the wave function of 
the electronic states at t = t* can be approximated as [2]: 



■*i(fi,f 2 ,...). (i) 



In this relation: i) ^(fi,^, ••■) is their wave function in the rest frame of the nucleus 
before the interaction; ii) r a is the coordinate vector of the a — th electron; iii) Qa = 
m e VA with m e electron mass. 

From cq. (1) the probability Pfi to reach the intrinsic state |^/ > starting from 
an initial state l^' > can be written as: 



fi = \< */ 



<*fl*'i >! 2 



< * / |e- i<5A -D- f, -|* j > 



(2) 



In the following, the case of the transition of an electron from one level to another 
and the case of its transition to the continuum (that is, the ionization of the atom) 
are separately analysed applying the reasonable approximation that multiple transi- 
tion/ionization processes can safely be neglected. 
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2.1 Calculation of the excitation probability and profile 

The total transition probability can be calculated from the single electron transition 
probability, Pj i , by exploiting the mean field approximation. In particular, we consider 
that the wave function of a single electron does not depend on the coordinates of the 
other ones. Thus, the probability Pj i of the transition of a single electron from an 
initial state \tpi > to a final bound state \ipf > is simply: 



P 1 



(3) 



where f is the coordinate vector of the considered electron. 

The term e~ t( ^ AT can safely be expanded in series by applying the dipole approxi- 
mation. In fact, since Va is at maximum of the order of the impinging WIMP velocity 
(less than the Galaxy escape velocity: <~ 2 • 10 _3 c) and R a ~ 10~ 10 m, we can write 
e = QaRo, — i^jMeVfm ~ ^ Thus, considering that in practice: ipi(r) ~ for r > R a 



and Qat < Qa t i the amplitude in cq. (3) can be re- written as: 

< ^ f \ e -iQ A ^. >c ,< ^|(i _ 2 Q Ar -_ i(Q A ff + ...)|^ > 



(4) 



retaining only the leading orders. 

As first, the probability for a single electron to retain the same state can be written 

as: 



P 1 = 



< 4>i\e 



\lpi > 



= F 



<iPi\{l-iQA?-...M> ~l-0(e 2 )<l, 



(5) 



where Fi is the atomic Form Factor of the orbital described by the ipi wave function. It 
can be derived from the Rayleigh scattering database (RTAB) [17], where the atomic 
Form Factors for the levels of the various atoms are reported. Some examples of the 
used atomic Form Factors for the Sodium and Iodine atoms are reported in Fig. 1; 
for the shell Is the atomic Form Factors calculated by means of the hydrogenic wave 
functions are also reported for comparison, showing a good agreement. 
On the other hand, for the transitions (/ ^ i) one gets: 



<ip f \(l-iQ A r- > ~ <ip f \Q A r\A> ~ 0(e 2 ) < 1 



(6) 



because of the states orthogonality. 

If we define 9 as the angle between the vectors Qa and r/i 
transition probability can be re-written as: 

P} i ^cos\6)Q\\r fi \ 2 . 



-< ipf\r\4>i >, this 



(7) 



Since the atoms are unpolarized, cos 2 (8) can be replaced by its average: 1/3, that is 



Pfi * \Q\ \rji 



(8) 
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Figure 1: Atomic Form Factors from [17] averaged over the shell for various principal 
quantum numbers of Sodium (left) and Iodine (right) atoms. For the shell Is the Form 
Factors calculated by means of the hydrogenic wave functions (dashed lines) are also 
reported for comparison. 

The Tfi matrix elements are connected with the standard dimensionless oscillator 
strengths, //j, by means of (as usual, here: % = 1): 

/,* = ^pM 2 (9) 

where uiu is the energy of the transition and, in conclusion, the transition probability 
can be expressed as: 



p i „ <pJfi_ (1Q) 
2m e ujfi 



The oscillator strengths for several atoms can also be derived from the Rayleigh 
scattering database [17]. Some examples of averaged oscillator strengths, calculated 
on this basis for various atoms, are reported in Table 1 together with some other 
evaluations available in literature. 



Let us now consider the general case. The probability of a single electron in the 
state i to have a transition to a free bound level (excitation) can be derived by summing 
over all the available possible states: 

nUi- E i<v^iqa^>i 2 , (ii) 

jGfree bound states 

neglecting the probability that other electrons fill the j-th state in the meantime. 
After an excitation, the atom returns to the ground state emitting X-rays and/or 
Auger electrons. The total energy released by the process is equal to the transition 
energy and it depends on the energy level of the final state. Considering that all the 
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Table 1: Some examples of the averaged oscillator strengths - calculated from the 
Rayleigh scattering database [17] - for various atoms and transitions, together with 
some other evaluations available in literature. 



Atom 


Transition 


//i from [17] 


ffi from literature 


H 


Is -> 2p 


0.416 


0.416 [18] 


H 


Is — > 3p 


0.079 


0.079 [18] 


H 


Is — » 4p 


0.029 


0.029 [18] 


Na 


3s -> 2p 


-0.042 


-0.043 [19] 


Na 


3s -> 3p 


0.968 


0.983 [19] 


Li 


2s -> 2p 


0.762 


0.753 [19] 


K 


4s — » 4p 


1.045 


1.02 [19] 



free bound levels are within at maximum few eV, we can safely neglect the spread of 
the excitation energies. Thus, the excitation profile of the z-th electron can simply be 
assumed to be a Dirac delta function centered around the average excitation energy, 
(E)i; this is generally less than few eV below the ionization threshold. Therefore, 
considering the excitations of the single i-th electron, the differential distribution of 
the electromagnetic part of the detected energy 2 , E em , for a given energy, Eq, provided 
by the WIMP to the nucleus at time t* , is given by: 



^A ex citation,-. 

Tie 



(E em \E ) = Pl und ME em - (E)i) . (12) 



2.2 Calculation of the ionization probability and profile 

The probability of a single electron in the state i to have a transition to the continuum 
(ionizing the atom), P} on j, can be derived as the difference between the unity and the 
probability of all the other possible processes. They are the following: i) the electron 
has no transition (P^ = F 2 ); ii) the electron inter-exchanges the states with another 
electron {P\ xc j); in) the electron excites to a free bound level {Pl ound Hence: 

pi _ 1 _ pi _ pi _ pi flO\ 

ion,i ii exc,i bound,i \ U J 

Taking into account the Pauli exclusion principle, the transition to an occupied 
state j is possible only if the j-th electron goes to another state. Thus, the exchange 
probability can be estimated by the expression: 

Plci* E \<^\Q A ^i>\ 2 {l-F]) (14) 

jfEoccupied bound states 

having neglected the probability that other electrons fill the^'-th state in the meantime. 

Unlike the excitation, the differential distribution, d i °^p t »°"' 1 ; considering the 
ionization due to the single i-th electron gives a continuum spectrum: 



ionization .i 



[E em \E ) = P^ MoniEem) , (15) 



dE e7n 

where the ionization profile D\ on {E em ) is evaluated in the following (see later). 



2 It is the sum of all the energies of the X-rays and Auger electrons, while in case of the ionization 
process - see later - also the energy of the electron escaping from the atom contributes. 
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Table 2: Estimated probabilities of excitation and ionization by Migdal effect dur- 
ing a 10 keV kinetic energy Na recoil or during a 33.3 kcV kinetic energy I recoil. 
Although the calculated probabilities are quite small, the unquenched nature of the 
electromagnetic contribution, the behaviour of the energy distribution for nuclear re- 
coils induced by WIMP-nuclcus elastic scatterings, etc. can give an appreciable impact 
at low WIMP masses; see later. 



Atom 


Shell 


Pu 


P 1 ■ 

exc.i 


P 1 

bound, i 


P 1 - 

ion, i 


< E > t (keV) 




Is 


0.99985 


1.7-10-'' 


7.8 • 10-'' 


1.5- 10" 4 


1.062 


Na 


2s 


0.99595 


1.2 • 10~ 5 


2.7 • 10" 5 


4.0 • 10~ 3 


0.062 




2p 


0.99595 


4.1 • 10" 6 


1.5 • 10" 4 


3.9- 10" 3 


0.033 




Is 


0.99999 


2.9- 10" iu 


2.8 • 10- 1U 


3.3- 10" b 


33.166 




2s 


0.99996 


8.4- 10~ 9 


7.3 • 10" 9 


4.5 • 10~ 5 


5.161 




2p 


0.99996 


9.8 • 10~ 9 


8.9 • lO- 9 


4.5 • 10~ 5 


4.687 




3s 


0.99966 


1.9 • 10~ 7 


1.1 • 10~ 7 


3.4- 10~ 4 


1.054 




3 P 


0.99966 


1.8 • 10~ 7 


1.5 • 10~ 7 


3.4- 10~ 4 


0.894 


I 


3d 


0.99966 


6.6 • 10~ 8 


3.5 • 10~ 8 


3.4- 10~ 4 


0.635 




4s 


0.99730 


5.5 • 10~ 6 


1.3 • 10~ 6 


2.7- 10~ 3 


0.190 




4p 


0.99730 


6.8 • 10~ 6 


3.8 • 10~ 6 


2.7- 10~ 3 


0.138 




4d 


0.99730 


6.5 • 10~ 6 


1.2 • 10" 5 


2.7- 10~ 3 


0.058 




5s 


0.98102 


2.5 • 10~ 4 


2.5 • 10" 5 


1.9- 10~ 2 


0.019 




5 P 


0.98102 


8.5 • 10~ 5 


7.3 • 10~ 3 


1.2- 10"2 


0.009 



In Table 2 some numerical estimates of the probabilities of the Migdal effect, dis- 
cussed in the previous and in the present sections, are given. It is worth to note that, 
although these probabilities are quite small, the electromagnetic unquenched nature of 
the contribution, the behaviour of the energy distribution for nuclear recoils induced 
by WIMP-nucleus clastic scatterings, etc. have some impact at low WIMP masses (see 
later). 

Considering the free electron approximation (that is, neglecting the interaction of 
the escaping electron with the nucleus), the final state of the escaping electron can be 

described as a plane wave [20]: ipp f ~ , normalized to a volume V. The final state 

, ., . Vd 3 p, 
density is: p f = -^f. 

Thus, the ionization probability, for an escaping electron with momentum pf, can 
be written as: 



< i>p f \e 



\tpi > 



Pf 



I 



-tpfl 



Vd 3 p 



f 



-i(QA+Pf)r 



-ipi(r)d 3 r 



d% = \Mp)\ 2 d 3 p f 



(2tt) 3 ~ 
Pi(p)d 3 p f , (16) 



where: i) p = pf + Qa] ii) <pi(p) is the wave function of the initial state in the 
momentum space; iii) pi (p) is the probability density to find a momentum p when the 
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electron is bound with wave function ipi. The ionization profile, p%{p), is normalized 
to 1, as it can be easily demonstrated: / 4>*(p)(j)i(p)d 3 pf — (<pi\<pi) = 1. 

Finally, for practical purposes Qa <))/; in fact, since the maximum energy pro- 
vided by the WIMP with mass, m^, and velocity, v, to a nucleus with mass, itla, 
is £ ,maz = \m w v 2 ^f^fp , the maximum velocity is: V A , ma x = J™+ mA v < 2v - 
Therefore, Qa <C 3 keV and can be safely neglected with respect to pf when the 
electron energy is larger than 10 eV. 

The energy distribution of the escaping electron can be obtained from eq. (16) by 

p 2 

introducing the variable Ej = (escaping energy of the electron) and by integrating 
over the solid angle. The angular integration can be easily performed by considering an 
average ionization profile for all the Z shell electrons in a full shell, where p s heii(Pf) — 
z \ n J2 s heii Pi(Pf) i s isotropic. Thus, we can define the ionization profile of a given 
shell as a function of the escaping energy of the electron according to: 

Dtf(Ef) ■ dE f ~ 4irm e ■ ^h^E s ■ p sheU ( Pf ) ■ dE f , (17) 

Moreover, in the real cases, the ionization is energetically allowed only for Ef, < 
Ef < Eq, with Eb binding energy of the electron. 




10 10 z 

p (keV) 



Figure 2: Example of pi s (p) and P2s{p) for Na atom as extracted from the Compton 
profile tables [21] by means of the procedure given in the text (stars). The solid 
lines - shown for comparison - are the p\ s (p) and p2 S (p) analytically calculated using 
hydrogenic wave functions and accounting for the screening effect of the inner electrons 
(Z e ff = 10 and Z e ff = 8 respectively). 



In order to obtain the functions p s heii(p), one can use the tables of the Compton 
profile Ji(p z ) — J J pi(p)dp x dp y available in literature [21] for the various shells of 
atoms 3 . As examples, pi s (p) and P2s(p) for Na atom are reported in Fig. 2 as 

3 In fact, for isotropic distributions (e.g. for s shell or for full shells) the following relation holds [22, 
23]: J she ii{pz) = J J Pshell(p)dpxdp y = 2n ppshell{p)dp and, hence: p shel i(p) = dJ shMv) _ 

Thus, the p s 

hellip) values of interest can be evaluated by means of the J s hell{p) tables given in ref. 

[21]. 
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extracted from the Compton profile tables [21]. They show a rather good agreement 
with the values calculated by using the hydrogenic wave functions and accounting for 
the screening effect of the inner electrons; these latter ones are also shown in the figure 
for comparison. 



3 Calculation of the expected counting rate 

For each given energy E , provided by the WIMP to the nucleus at time t* , the 
differential distribution, J i E Wi , considering the transitions (ionizations/excitations) of 
the single i-th electron, is given by the sum of the two contributions of cq. (12) and eq. 
(15). Moreover, the probability density distribution to have an electromagnetic release, 
E em , for a given E value (considering all the contributions from all the electrons in 
the atom) can be written as: 



^i(E em \E ) ~ P Eo S(E em ) +Y J ^^(E em \E,) + ... , (1<S) 



where the term P Eo is the probability for the whole atom to remain unchanged: P E ~ 
Eli P i,E > and p Ie q = 1 - P ion,i ~ P bound,i is the probability that a single electron has no 
transition either to free bound level or to the continuum. In eq. (18) the expansion can 
be stopped at the shown level because the probability of multiple excitation/ionization 
is negligible. 

Finally, let us rewrite eq. (18) as a function of the detected energy, Edet, which 
is given by the sum of the recoil detected energy, Ef r , and of the electromagnetic 
component (the relation + E em — Eq holds): 

<lN -(E det \E Q ) = ^(E em \E ) ■^ = T ^—- ^(E em \E ) , (19) 



dEdet dE em dEdet 1 — Q_a dE e 

where qA is the nuclear recoil quenching factor for the considered nucleus in the given 
detector at the considered energy 

In particular, let us now point out the case of recoils induced by WIMP-nucleus 
elastic scatterings under the usual hypothesis that just one component of the dark 
halo can produce elastic scatterings on nuclei. 

For every target specie A, the expected energy distribution including the Migdal 

dR {M) 

effect, dE A d - (E d et), requires the Eq differential distribution produced in the WIMP- 
nucleus elastic scattering, given in squared brackets: 



dE de 



(E det )= I -^(E det \E ) N T ^ £r(v,Eo)vf(v)dv 
J dEdet m w Jv min (E ) dE 



dE .(20) 



There: i) Nt is the number of target nuclei of A specie; ii) pw — (,Po, where p is the 
local halo density and £ < 1 is the fractional amount of local WIMP density; iii) f(v) 
is the WIMP velocity (v) distribution in the Earth frame: iv) v m i n = / " 1 ^- E " (raw A 

V WA 

is the reduced mass of the WIMP-nucleus system) ; v) v max is the maximal WIMP ve- 
locity in the halo evaluated in the Earth frame; vi) j§^(v, E ) = ( ) + ( ) , 

V / SI V / SD 
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with ( -jj^ ) spin independent (SI) contribution and ( -jjjr ) spin dependent (SD) 
contribution. 

Finally, the expected differential counting rate as a function of the detected energy, 
Edet, for a real multiple-nuclei detector (as e.g. the Nal(Tl)) can be easily derived by 
summing the eq. (20) over the nuclei species and taking into account the detector 
energy resolution: 

aet J A=Na,I 

The G(Ed et , E') kernel generally has a gaussian behaviour. 

Obviously the expected differential counting rate has to be evaluated in given 
astrophysical, nuclear and particle physics scenarios, also requiring assumptions on all 
the parameters needed in the calculations and the proper consideration of the related 
uncertainties (for some discussions see e.g. [4, 5, 6]). 
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Figure 3: Examples of shapes of expected energy distributions from WIMP-nucleus 
elastic scatterings in the Nal(Tl) detectors of DAMA/Nal with (continuous line) and 
without (dashed line) including the Migdal effect, for the model framework given in 
the text. The effect of the inclusion of this existing physical effect is evident. The 
vertical lines indicate the DAMA/Nal energy threshold. 



For clarity, Fig. 3 shows just few examples of shapes of expected energy distribu- 
tions with and without accounting for the Migdal effect. For this template purpose - 
accounting also for the experimental features of the detectors [7, 8, 4, 5] - we just adopt 
the following assumptions among all the possibilities: i) WIMP with dominant Spin 
Independent coupling and with nuclear cross sections oc A 2 ; ii) non- rotating Evans' 



logarithmic galactic halo model with core radius R c = 5 kpc, local velocity vo = 170 
km/s and po = 0.42 GeV cm~ 3 (Bl halo model in rcf. [4, 5]); iii) form factors and 
quenching factors of 23 Na and 127 I as in case C of ref. [4]. The used normalizations 
assure the same vertical scale in the shown plots. It is clear the fraction of events 
at very low WIMP masses of electromagnetic nature. Note that other choices of the 
model framework do not change the substance of the results. 

4 Some examples 

The proper accounting of the Migdal effect in corollary quests for WIMPs as DM 
candidate particles can be investigated by exploiting the expected energy distribution, 
derived above, to some of the previous analyses on the DAMA/Nal annual modulation 
data in terms of WIMP-nucleus elastic scattering. For this purpose, the same scaling 
laws and astrophysical, nuclear and particles physics frameworks of refs. [4, 5] are 
adopted, while - for simplicity to point out just the impact of the Migdal effect - the 
SagDEG contribution to the galactic halo, whose effect we discussed in ref. [6], will 
not be included here. 

The results for each kind of interaction are presented in terms of allowed vol- 
umes/regions, obtained as superposition of the configurations corresponding to likeli- 
hood function values distant more than 4tr from the null hypothesis (absence of mod- 
ulation) in each one of the several (but still a very limited number) of the considered 
model frameworks. This allows us to account - at some extent - for at least some of 
the existing theoretical and experimental uncertainties (see e.g. in ref. [4, 5, 9, 6] and 
in the related astrophysics, nuclear and particle physics literature). 

Since the 23 Na and 127 I are fully sensitive both to SI and to SD interactions, the 
most general case is defined in a four-dimensional space (mw, £,&si, fi&SD, 0), where: 
i) a si is the point-like SI WIMP-nucleon cross section and <jsd is the point-like SD 
WIMP-nucleon cross section, according to the definitions and scaling laws considered 
in ref. [4] ; ii) tg6 is the ratio between the effective coupling strengths to neutron and 
proton for the SD couplings (9 can vary between and n) [4] . The subcase of purely SI 
coupled WIMPs is shown in Fig. 4, while in Fig. 5 just two slices of the 3-dimensional 
allowed volume (mw, £csd, 0) for the purely SD case are given as an example. 

It is worth to note that the accounting for the electromagnetic aspects of the interac- 
tions provides in the considered scenarios at the given C.L. additional volumes/regions 
not topologically connected with the remaining allowed parts. This depends on the 
behaviour of the expected energy distributions at low masses (where the Migdal con- 
tribution is appreciable) with respect to that at higher masses, where recoils dominate. 

Finally, in the general case of mixed SI&SD coupling one gets a 4-dimensional 
allowed volume (£05/ , ^<jsd, mw, &)\ new allowed volume at the given C.L. is present 
in the GeV region. Fig. 6 shows few slices of such a volume as examples. 

Note that general comments, extensions and comparisons already discussed in ref. 
[4, 5, 9, 6] still hold. 

Finally, just for completeness, we remind that GeV mass DM particles have been 
proposed in ref. [24, 25, 26] in order to offer a mechanism able to account for the Baryon 
Asymmetry in the Universe and to naturally explain why Qdm ~ 5Q;,. Moreover, in 
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10 

m w (GeV) 



Figure 4: Region allowed in the (£crsi, mw) plane in the considered model frameworks 
for pure SI coupling; see text. The hatched region appears when accounting for the 
Migdal effect. Inclusion of other contributions and/or of other uncertainties on param- 
eters and models, such as e.g. the SagDEG contribution [6] or more favourable form 
factors, would further extend the region and increases the sets of the best fit values. 
For completeness and more see also [4, 5, 9, 6]. 




m w (GeV) m w (GeV) 

Figure 5: Two slices of the 3-dimensional allowed volume (£<jsd, mw, 0) m the con- 
sidered model frameworks for pure SD coupling; see text. The hatched regions appear 
when accounting for the Migdal effect. Analogous remarks as those in the caption of 
Fig. 4 hold. 

ref. [24] it was shown that a GeV mass DM candidate would potentially solve the 
discrepancies between observations and KCDM model on the small scale structure 
of the Universe. Finally, among the GeV mass WIMP candidates we remind: i) the 
H dibaryon, already predicted within the Standard Model of particle Physics [26]; 
ii) the Darkon, a real scalar field in an extended Standard Model [27]; iii) the light 
photino early proposed in models of low-energy supersymmetry [28] ; iv) the very light 
neutralino in Next-to-MSSM model [29]; v) the scalar GeV mass DM candidates of 
ref. [30]; vi) the mirror Deuterium in frameworks where mirror matter interactions 
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Figure 6: Examples of slices of the 4-dimensional allowed volume (£,<jsi, £,&SD, mw, 6) 
in the considered model frameworks; see text. The hatched slice appears when ac- 
counting for the Migdal effect. Analogous remarks as those in the caption of Fig. 4 
hold. 



[31] with ordinary matter are dominated by very heavy particles. 



5 Conclusions 



In this paper the ionization and the excitation of bound atomic electrons induced by 
the perturbation of the recoiling nucleus after a WIMP elastic scattering have been 
discussed. This effect has so far usually been neglected in the field. The needed 
theoretical arguments have been developed and the related impact in corollary quests 
for the candidate particle has been shown, as example, for some simplified scenarios. 
Obviously, many other arguments can be addressed as well both on DM candidate 
particles and on astrophysical, nuclear and particle physics aspects; for more see [4, 5, 
9, 6] and in literature. 
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